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Here,  we  gifce  the  final  report  for  AFOSR  Contract  85-012&-CfTindamental 
Aspects  of  the  Structure  of  Supersonic  Turbulent  Boundary  Layers, "^monitored  by 
Dr.  J.  McMichael.  This  three-year  contract  originally  had  two  tasks:  to 
investigate  the  structure  of  flat  plate  supersonic  turbulent  boundary  layers  ^ 

(Task  A),  and  to  investigate  longitudinal  curvature  effects  in  turbulent  boundary 
layers  (Task  B).  As  a  result  of  the  work  performed  under  these  task  headings,  we \ 
identified  a  need  to  study  the  structure  of  simple,  wall-bounded  vortex  loops.  \ 
This  investigation  (Task  C)  commenced  in  the  third  year,  and  we  feel  that  it  has  \ 
made  a  major  contribution  to  our  understanding  of  the  ensemble-averaged  \ 

structure  of  low  and  high  speed  turbulent  boundary  layers.  The  progress  \ 

reported  here  is  that  achieved  in  completing  these  three  tasks. since  the  \ 


beginning  of  the  Contract  in  January  1985. 
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INTRODUCTION 


Here,  we  give  the  final  report  for  AFOSR  Contract  85-0126  "Fundamental 
Aspects  of  the  Structure  of  Supersonic  Turbulent  Boundary  Layers,”  monitored  by 
Dr.  J.  McMichael.  This  three-year  contract  originally  had  two  tasks:  to 
investigate  the  structure  of  flat  plate  supersonic  turbulent  boundary  layers 
(Task  A),  and  to  investigate  longitudinal  curvature  effects  in  turbulent  boundary 
layers  (Task  B).  As  a  result  of  the  work  performed  under  these  task  headings,  we 
identified  a  need  to  study  the  structure  of  simple,  wall-bounded  vortex  loops. 
This  investigation  (Task  C)  commenced  in  the  third  year,  and  we  feel  that  it  has 
made  a  major  contribution  to  our  understanding  of  the  ensemble-averaged 
structure  of  low  and  high  speed  turbulent  boundary  layers.  The  progress 
reported  here  is  that  achieved  in  completing  these  three  tasks  since  the 
beginning  of  the  Contract  in  January  1985. 
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TASK  A:  THE  STRUCTURE  OF  SUPERSONIC.  TURBULENT  BOUNDARY  LAYERS 


The  principal  aim  of  this  task  was  to  develop  a  physical  model  for 
supersonic  turbulent  boundary  layer  structure  which  could  be  used  as  the  basis 
for  improved  calculations  of  turbulent  flow.  The  effort  was  specifically  directed 
towards  a  better  understanding  of  the  large-scale  structures,  that  is,  the 
energy-containing  motions  in  the  fully  turbulent  region,  under  conditions  of  zero 
and  adverse  pressure  gradients. 

The  presence  of  "finger-like"  structures  in  a  zero  pressure  gradient 
boundary  layer  was  established  by  Spina  &  Smits  (1987).  These  structures  are 
inclined  at  about  45°  to  the  flow  direction,  they  appear  to  have  a  limited 
spanwise  and  streamwise  extent,  and  they  are  convected  downstream  for  consider¬ 
able  distances  without  losing  their  identity.  These  conclusions  were  based  on 
measurements  taken  with  multiple  normal  hot-wire  probes,  multiple  wall-pressure 
transducers,  and  high-contrast  schlieren  photographs.  The  analysis  used 
space-time  correlations  and  conditional  averaging  using  the  single-point  VITA 
technique. 

Some  interesting  new  results  have  been  obtained  since  then. 

(a)  Eddy-angle  measurements 

Measurements  of  the  structure  angle  (derived  from  the  space-time  correlation 
of  the  signals  from  two  parallel  wires)  have  been  obtained  for  an  adverse 
pressure  gradient  supersonic  layer  and  a  zero  pressure  gradient  subsonic  layer. 

In  the  supersonic  case,  the  initial  Reynolds  number  based  on  momentum  thickness 
was  about  82,000,  the  static  pressure  rose  by  a  factor  of  1.9  in  a  distance  of 
approximately  10  6o  (do  is  the  initial  boundary  layer  thickness),  and  the  Mach 
number  decreased  from  2.85  to  2.10  in  the  same  distance  (Fernando  and  Smits, 


1987).  The  results  are  shown  in  Fig.  1  for  two  different  wire  separations.  The 
results  of  the  smaller  separation  wires  (t/6o  =  0.09)  show  that  the  incoming 
boundary  layer  structure  angle  is  around  45°  over  most  of  the  layer  with  a 
decrease  in  angle  near  the  wall  and  an  increase  in  angle  near  the  outer  edge  of 
the  layer.  The  structure  angle  can  be  seen  to  increase  by  approximately  5° 
through  the  interaction.  Much  the  same  results  were  obtained  with  the  larger 
separation  wires  (t /6c  =  0.18)  although  the  structure  angles  are  consistently 
higher  by  around  5° .  These  angles  were  calculated  assuming  a  convection  velocity 
equal  to  the  mean  velocity  half-way  between  the  wires.  The  data  have  recently 
been  re-analyzed  assuming  a  fixed  (invariant  with  wall  distance)  structure 
convection  velocity  of  0.8  times  the  free-stream  velocity.  The  differences  in 
structure  angles  between  the  incoming  and  the  adverse  pressure  gradient  boundary 
layers  were  found  to  be  slightly  smaller  and  of  the  order  of  the  experimental 
uncertainty. 

For  the  subsonic  flow,  structure  angles  have  also  been  derived  from  the 
space-time  correlation  of  signals  from  two  normal  wires  (i.e.  the  technique  of 
Spina  &  Smits,  1987).  The  results  were  obtained  for  a  flat  plate,  zero  pressure 
gradient  layer,  and  a  boundary  layer  recovering  from  strong  convex  curvature 
(Alving  &  Smits  1986;  Alving  et.al.  1987).  The  results  for  the  zero  pressure 
gradient  boundary  layer  with  R^  =  5700  show  that  the  angle  lies  between  40  and  50 
degrees  for  y/6o  <  0.5  (see  Fig.  1).  However,  the  angle  increases  further  from 
the  wall  and  reaches  a  value  around  60°  near  the  edge  of  the  layer.  There 
appears  to  be  a  consistent  increase  in  9  following  the  convex  curvature  which  may 
be  related  to  the  new  equilibrium  state  observed  in  the  Reynolds  stress  profiles. 
This  may  imply  that  several  different  types  of  events  may  be  contributing  to  the 
space-time  correlation,  and  that  some  type  of  conditional  sampling  will  probably 
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be  necessary  to  sort  out  the  influence  of  the  various  types  of  events.  Work  is 
continuing  in  this  area. 


Two  additional  results  have  been  obtained  with  respect  to  the  eddy  angle  in 
supersonic  flow.  First,  by  using  VITA  to  detect  significant  events  on  one  wire, 
and  examining  the  simultaneous  signal  from  the  other  wire,  the  variance  of  the 
eddy  angle  can  be  determined.  Preliminary  results  have  shown  that  the  two 
signals  typically  have  a  correlation  coefficient  greater  than  0.95,  and  that  the 
variance  is  of  the  order  of  15° -  Second,  one  wire  was  placed  as  close  as 
possible  to  the  wall,  at  y*  e  60,  and  another  wire  was  traversed  through  the 
layer.  This  method  of  measuring  eddy  angles  has  been  tried  by  others,  and  gives 
very  low  angles  throughout  the  layer,  of  order  15°  to  30°  (Brown  &  Thomas  1977, 
Robinson  1986).  Our  measurements  gave  similar  results.  As  pointed  out  by  Spina  & 
Smits  (1987),  however,  the  variation  of  convection  velocity  with  distance  from 
the  wall  makes  these  measurements  difficult  to  interpret,  and  they  appear  to  be 
of  limited  value.  Using  two  wires  placed  a  short  distance  apart  to  deduce 
structure  angles  is  physically  more  meaningful. 

(b)  Conditionally  Sampled  Results 

In  our  previous  work,  the  primary  basis  for  conditional  sampling  was 
the  VITA  technique,  and  it  was  applied  to  signals  from  a  normal  wire  which 
measures  mass-flow  fluctuations  (pu)’.  The  (pu)’  signal  was  converted  to  u’  using 
the  Strong  Reynolds  Analogy.  The  results  obtained  in  zero  pressure-gradient 
high-Reynolds  number  supersonic  flows  agreed  well  with  similar  results  from 
subsonic,  low  Reynolds  number  flows.  However,  it  is  not  at  all  clear  what 
physical  meaning  can  be  attributed  to  the  events  detected  by  this  technique.  The 
interpretation  is  limited  by  the  fact  that  VITA  is  a  single  point  criterion, 
applied  to  measurements  of  a  single  component  of  the  turbulent  fluctuations.  Our 
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current  emphasis,  therefore,  is  to  improve  the  detection  technique  to  the  point 
where  the  detected  events  can  be  given  a  more  complete  physical  interpretation. 

We  have  already  found  that  the  VITA  technique  can  give  a  distorted  picture 
of  the  flow  structure.  The  primary  technique  which  has  allowed  this  critical 
interpretation  of  VITA  is  the  development  of  a  reliable  and  accurate  crossed- 
wire  probe,  capable  of  giving  signals  proportional  to  instantaneous  values  of  u’, 
v’  (and  u’v’)  (see  Fernando,  Donovan  &  Smits,  1987).  The  shear  stress  signal  u’v’ 
is  highly  intermittent,  with  sharp  excursions  from  the  background  level  (see  Fig. 
2).  This  signal  is  extremely  suitable  for  conditional  sampling.  A  simple 
threshold  criterion  was  used  to  detect  the  events  which  make  the  greatest 
contribution  to  shear  stress  (for  this  case  the  threshold  =  2  u’ v’ ) .  These  events 
are  obviously  of  great  practical  interest.  Quadrant  analysis  was  used  to  classify 
the  u’v’  events  into  four  categories:  I  (u+,  v+),  II  (u-,  v+),  III  (u-,  v-),  and 
IV  (u+,  v-). 

The  results  were  very  revealing.  Typically,  the  ensemble-averaged  events 
were  strong,  and  rather  simple  in  shape  (see  Fig.  3).  In  particular,  it  should  be 
noted  that  the  u’  event  is  single-sided,  that  is,  it  is  either  positive  or 
negative.  In  contrast,  when  VITA  was  applied  to  the  u’  signal  (instead  of 
thresholding  on  u’v’),  the  event  was  doubled-sided,  as  shown  in  Fig.  4.  These 
results  suggest  that  VITA  detection  is  subject  to  ambiguity,  probably  caused  by 
the  superposition  of  two  types  of  events  which  can  only  be  separately  identified 
by  using  quadrant  analysis. 

To  explain  the  concept,  consider  two  individual  u  events  detected  by  the 
shear  stress  criteria,  one  lying  in  quadrant  II  and  the  other  in  quadrant  IV  of 
the  u-v  plane.  Both  are  single  sided  events.  It  is  possible  that  VITA  detects 
both  these  events  as  positive  by  checking  the  slope  at  their  centers.  The 


ensemble  average  of  these  two  events  would  be  a  two  sided  event.  Figure  5 
clarifies  the  concept.  Depending  on  exactly  where  the  VITA  technique  checks  the 
slope,  thse  events  could  have  been  picked  up  a  negative  events,  giving  a  two 
sided  negative  event.  This  may  explain  why  positive  and  negative  VITA  events  are 
nearly  mirror  reflections  of  each  other,  about  the  velocity  axis. 

The  quadrant  analysis  can  be  used  to  reveal  another  interesting  result.  The 
ensemble-averaged  events  in  each  quadrant  (as  shown,  for  example,  in  Fig.  3)  plot 
trajectories  in  the  u-v  plane.  Some  examples  of  these  loci  are  shown  in  Fig.  6. 
These  trajectories  will  be  compared  to  the  results  of  the  inviscid  vortex 
skeleton  model  calculations  described  in  Task  C  Progress.  The  probability 
density  of  the  shear  stress  in  the  u-v  plane  for  the  upstream  boundary  layer  is 
plotted  for  various  y/6o  in  Fig.  7.  There  is  a  definite  change  in  the  character 
of  the  probability  density  as  y/60  increases.  In  the  lower  regions  of  the 
boundary  layer  all  the  equi-probability  density  contours  are  shaped  more  or  less 
like  ellipses,  centered  around  the  origin  of  the  u-v  plane,  with  the  major  axis 
lying  in  quadrants  II  and  IV  at  an  angle  to  the  u  axis.  As  y/6o  increases,  the 
major  axis  rotates  and  aligns  itself  with  the  u  axis.  The  centers  of  the 
equi-probability  density  contours  show  two  different  effects:  the  higher 
equi-probability  contour  centers  move  into  quadrant  IV,  and  the  remaining  contour 
centers  move  into  quadrant  II.  This  behavior  of  the  equiprobability  contours 
with  increasing  y/6o  does  not  change  within  the  interaction.  The  contribution  to 
the  shear  stress  from  each  quadrant  and  the  fractional  residence  time  in  each 
quadrant  is  shown  in  Figs.  8  and  9  respectively,  for  various  x  positions.  It  is 
seen  that  the  trend  in  the  fractional  contribution  to  shear  stress  with  y/<50  is 
the  same  at  all  stations.  Lower  in  the  boundary  layer,  contributions  from 
quadrants  II  and  IV  are  approximately  the  same.  Higher  in  the  boundary  layer  the 
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contribution  from  quadrant  II  is  larger,  even  though  the  fractional  residence 
time  in  quadrant  II  is  smaller.  This  implies  that  the  contributions  to  the  shear 
stress  from  quadrant  II  consists  of  more  intense  shorter  duration  peaks  in  u’v*. 
This  same  conclusion  can  be  drawn  from  the  "Hole  Analysis"  suggested  by  Lu  and 
Willmarth  (1970),  and  by  examining  the  probability  density  of  the  time  spent 
above  a  certain  threshold  shear  stress,  in  quadrants  II  and  IV.  At  all  y/S0 
there  is  a  finite  negative  contribution  to  the  shear  stress  from  quadrants  I  and 
III.  Hence,  the  contributions  from  quadrants  II  and  IV  add  up  to  more  than  the 
total  shear  stress.  As  can  be  seen,  in  the  outer  part  of  the  boundary  layer,  the 
contribution  to  the  shear  stress  from  quadrant  II  alone  can  be  greater  than  the 
time-averaged  shear  stress. 

c)  Conditionally-Sampled  Flow  Visualization 

We  have  continued  to  improve  our  flow  visualization  techniques.  In 
particular,  we  are  continuing  to  refine  a  method  of  conditionally  sampling 
schlieren  images  of  boundary  layer  structure  in  supersonic  flow.  In  this 
technique,  a  hot-wire  probe  is  used  to  detect  the  presence  of  strong,  large- 
scale  motions  by  using  a  real-time  analog  of  the  VITA  detection  method.  Upon 
detection,  a  light  source  is  flashed  to  record  a  microsecond  exposure  schlieren 
image  on  video  tape.  The  hot-wire  output  is  recorded  simultaneously  on  the  same 
video  frame,  and  a  typical  image  is  shown  in  Fig.  10.  The  image  was  recorded 
using  a  CCD  camera  with  linear  gain  so  that  the  intensity  can  be  linearly  related 
to  density  gradient,  and  images  can  be  added  to  get  an  ensemble-averaged  picture. 
The  conditionally  sampled  nature  of  the  images  means  that  we  record  only  large 
events  which  are  in  the  plane  of  the  wire.  We  hope  this  method  will  remove  the 
spatial  integration  effect  of  schlieren  images.  Preliminary  results  were 
presented  by  Smith  and  Smits  (1986). 


We  have  recently  made  high  speed  (40K  frames/sec)  laser  schlieren  movies  of 
a  supersonic  boundary  layer  which  clearly  show  the  presence  of  strong  large  scale 
motions  convecting  with  the  layer(Smith  &  Smits  1988).  Figure  11  is  a  typical 
sequence  of  frames  taken  from  one  of  these  movies.  It  is  possible  to  "threshold" 
the  schlieren  images  so  that  only  the  very  strongest  structures  are  detected. 
Despite  the  spatial  integration  inherent  in  the  schlieren  technique  it  appears 
possible  to  produce  convincing  images  of  individual  structures  and  their 
evolution. 


In  our  current  work,  we  are  using  a  rake  of  8  hot-wires  spaced  in  the 
vertical  direction,  and  a  probe  consisting  of  6  hot-wires  arranged  in  the  cross - 
stream  direction,  to  characterize  the  large-scale  motions  in  more  detail.  The 
results  are  under  analysis. 


TASK  B:  LONGITUDINAL  CURVATURE  EFFECTS  IN  TURBULENT  BOUNDARY  LAYERS 


Here,  we  study  the  distortion  and  relaxation  of  an  initially  self-preserving 
zero-pressure  gradient  subsonic  turbulent  boundary  layer  as  it  passes  through  a 
region  of  convex  curvature  (see  Fig.  12).  The  study  serves  two  major  purposes: 
it  provides  us  with  measurements  in  a  subsonic  zero-pressure  gradient  layer  which 
can  be  used  for  direct  comparison  with  the  corresponding  measurements  in  the 
supersonic  case,  and  it  allows  us  to  study  the  behavior  of  structural  parameters 
subsequent  to  the  severe  distortion  imposed  by  the  convex  curvature. 
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Previously,  we  were  concerned  with  achieving  an  acceptable  level  of  two- 
dimensionality  in  the  relaxing  boundary  layer,  while  at  the  same  time  minimizing 
the  pressure  gradients  in  the  bend.  We  believe  these  aims  were  achieved,  and  we 
are  satisfied  with  the  overall  flow  behavior. 

A  full  set  of  Reynolds  stress  measurements  has  been  completed,  and  the 
results  were  presented  at  the  25th  Aerospace  Sciences  Meeting  in  Reno  (Alving, 
Watmuff  and  Smits  1987).  The  results  are  shown  in  Fig.  13.  The  collapse  of  the 
shear  stress  in  the  outer  layer  at  the  exit  of  the  bend,  and  the  subsequent 
recovery  behavior  is  clearly  shown.  At  station  11,  about  1206c.  downstream  of  the 
end  of  curvature,  the  recovery  of  the  shear  stress  profile  is  obviously  far  from 
complete. 

The  results  were  obtained  using  dynamic  calibration.  The  advantage  of  this 
method  is  that  the  sensitivity  to  both  velocity  components  can  be  determined 
directly  without  the  need  to  measure  wire  angles  or  assume  heat  transfer  laws.  In 
the  method  described  by  Watmuff  et  al.  (1983),  the  out-put  signals  from  a 
crossed-wire  probe  were  combined  to  produce  signals  sensitive  only  to  the 
streamwise  and  the  normal  velocity  components  respectively.  This  "matching" 
procedure  was  performed  by  adjusting  analog  computer  circuits  while  the  probe  was 
subjected  to  velocity  perturbations  in  the  free  stream  by  a  mechanical  shaker. 

The  probe  was  then  oscillated  once  more  in  both  directions  and  phase-averaged 
values  of  the  matched  voltages  were  used  to  determine  the  sensitivities  at  a 
number  of  free-stream  velocities.  The  final  calibration  was  obtained  by 
analytically  integrating  polynomial  expressions  curved-fitted  to  the 
sensitivities.  Here,  we  have  developed  a  new  technique  where  the  matching  can  be 
calculated  from  the  phase-averaged  values  of  the  anemometer  signals.  The 


sensitivities  of  the  matched  signals  can  be  then  be  calculated  from  the  original 
data  to  produce  the  final  calibration  without  the  need  to  impose  further  velocity 
perturbations.  The  method  is  considerably  less  time  consuming  than  that  used  by 
Watmuff  et  al.  since  the  probe  does  not  have  to  be  oscillated  both  directions 
twice  (once  for  matching  and  then  for  final  calibration)  and  it  also  removes  any 
operator  bias  in  the  matching  of  the  wires. 

Our  more  recent  focus  for  Task  B  has  meant  that  instead  of  concentrating  on 
Reynolds  stress  measurements  for  a  range  of  curvature  parameters,  we  have  carried 
out  measurements  of  structural  parameters.  This  approach  allows  a  more 
fundamental  understanding  of  curvature  effects,  and  provide  a  direct  comparison 
with  the  supersonic  boundary  layer  results.  The  first  of  these  comparisons  was 
made  by  measuring  the  structure  angle  in  the  self-preserving  boundary  layer 
upstream  of  curvature.  The  results  were  shown  in  Fig.  1,  and  they  were  discussed 
under  Task  A.  Measurements  have  also  been  made  using  the  rake  of  8  hot  wires 
described  under  Task  A.  The  results  are  currently  being  analyzed  and  will  be 
presented  by  Alving  and  Smits  (1988). 

TASK  C:  ANALYSIS  AND  CONTROL  OF  WALL-BOUNDED  VORTEX  LOOPS 

It  has  been  widely  suggested  that  the  structure  of  turbulent  boundary  layers 
can  be  described  in  terms  of  "hairpin"  or  "lambda"-shaped  vortex  loops.  This 
physical  model  seems  to  describe  equally  well  the  dynamics  of  small-  and  large- 
scale  motions  (see,  for  example,  Head  &  Bandyophadyay ,  1981,  and  Smith,  1984). 

The  interpretation  of  our  conditionally-sampled  measurements  has  been  strongly 
influenced  by  the  hairpin  loop  model,  and  to  aid  this  interpretation  further  we 
have  begun  to  investigate  the  behavior  of  simple  wall-bounded  vortex  loops. 


At  first,  we  tried  generating  these  loops  using  the  vortices  shed  from  a 
hemisphere  placed  on  the  wall,  similar  to  that  used  by  Perry  et  al.  (1982)  and 
Acarlar  &  Smith  (1984).  This  method  proved  not  to  be  suitable  for  our  purposes 
because  we  wished  to  control  the  frequency  of  generation.  By  controlling  their 
formation,  interactions  among  multiple  vortex  loops  can  be  studied.  The  method 
currently  in  use  employs  a  lambda-shaped  wire  which  periodically  forces  laminar, 
vortical  fluid  away  from  the  wall  (see  Fig.  14).  The  experiments  were  conducted 
in  a  large  low-speed  smoke  tunnel  so  that  combined  flow-visualization  and  hot¬ 
wire  measurements  were  possible. 


Task  C  Progress 

The  nature  of  the  single  loops  is  illustrated  in  Figs.  15  and  16.  Figure  15 
shows  a  plan  view  of  the  loop,  indicating  the  characteristic  splaying  of  the  legs 
of  the  vortex,  and  the  omega-shaped  head.  Figure  16  shows  a  side  view,  together 
with  signals  proportional  to  u’  and  v’,  as  measured  by  the  hot-wire  probe  visible 
on  the  left.  The  signals  are  similar  to  the  ensemble-averaged,  sec^rd-quadrant  u* 
and  v’  signals  shown  in  Fig.  3. 

Simple  analytical  models  have  been  developed  to  model  these  experimental 
vortex  loops  using  straight-line  vortex  filaments.  The  example  given  in  Fig.  17 
shows  a  simple  line  array  of  five  loops,  and  the  u’  and  v*  signals  corresponding 
to  two  positions  are  shown  in  Fig.  18.  The  similarities  to  the  ensemble-averaged, 
second-  and  fourth-quadrant  u’  and  v*  signals  shown  in  Fig.  3  are  obvious.  The 
loci  on  the  (u’  :  v’)  plane  for  different  observers  are  shown  in  Fig.  19,  and 
there  are  some  interesting  comparisons  possible  with  the  results  shown  in  Fig.  6. 

We  have  continued  to  use  the  vortex  loop  study  as  support  for  the 
interpretation  of  our  conditionally  sampled  measurements.  We  have  commenced  with 
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a  study  of  interacting  loops.  We  have  already  performed  some  preliminary  work  to 
examine  vortex  pairing  by  generating  two  different  strength  vortex  loops  down  the 
center  of  the  tunnel,  and  demonstrated  that  vortex  pairing  in  wall-bounded  flows 
is  an  unlikely  mechanism  for  the  growth  of  the  attached  eddies. 


COMPLEMENTARY  WORK 

Prof.  A.  E.  Perry  has  visited  on  two  occasions,  (four  weeks  in  January  1986 
and  two  weeks  in  January  1987).  His  suggestions  and  criticisms  have  proven  to  be 
very  useful  in  designing  our  experimental  work.  Our  current  cooperation  is 
centered  on  developing  attached  eddy  models  for  compressible  boundary  layers. 

In  June  1986,  we  invited  Stephen  K.  Robinson  from  NASA-Ames  to  visit  us  and 
discuss  supersonic  turbulent  boundary  layer  behavior.  Steve  Robinson  is  one  of 
the  very  few  people  who  have  approached  supersonic  boundary  layer  measurements 
with  a  view  to  understanding  the  structure  of  the  flow  (see,  for  example, 
Robinson,  1986).  The  visit  was  very  constructive,  and  we  followed  up  his  1986 
visit  with  another  in  1987.  At  that  time,  Robinson  brought  along  a  hot-wire  rake 
of  8  wires  which  we  implemented  in  a  zero  pressure  gradient  supersonic  flow  with 
his  active  cooperation.  The  results  are  being  analyzed. 

The  meetings  of  the  Workshop  on  High-Speed  Boundary  Layers  have  continued. 
Meetings  were  held  to  coincide  with  the  Atlanta  and  Reno  AIAA  meetings,  and  a 
special  workshop  on  shock-wave/boundary  layer  interactions  was  held  at  Penn-State 
in  July  1987. 

The  list  of  publications  produced  with  support  of  AFOSR  Grant  85-0126  is 
given  in  Appendix  A. 
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CONCLUSIONS 

The  work  performed  under  this  contract  has  succeeded  on  identifying  and 
describing  the  large  scale,  coherent  motions  in  a  Mach  3  turbulent  boundary 
layer.  Both  zero  and  adverse  pressure  gradient  flows  were  investigated.  A 

*4 

complementary  study  examined  the  relaxation  a  boundary  layer  downstream  of  a 

A 

short  region  of  convex  curvature,  in  terms  of  the  behavior  of  the  large-scale 
motions.  The  results  of  both  investigations  have  been  interpreted  in  terms  of 
vortex  loop  models  of  turbulent  boundary  layers,  and  have  provided  the 
preliminary  basis  for  the  development  of  predictive  models  incorporating  concepts 
based  on  coherent  motions. 
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Signals  of  u’v',  v’  and  u’  deduced  from  the  output  of  a 
crossed-wire  in  zero  pressure  gradient  supersonic  boundary  layer 
at  y/6  =  0.65. 
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Fig.  10.  Microsecond  schlieren  video  image  of  zero  pressure  gradient 

supersonic  boundary  layer.  The  hot-wire  probe  is  on  the  left,  and 
the  corresponding  signal  is  shown  near  the  bottom  of  the  picture. 
The  time  scale  corresponds  approximately  to  the  length  scale  of  the 
picture. 
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Fig.  11.  High-speed  laser  schlieren  movie  of  zero  pressure  gradient  supersonic 
boundary  layer.  The  images  havp  been  digitally  processed  to  enhance 

the  contrast .  Sequence  of  poeitivt  event!  TV  Art!  from*  i»  In  the  epper  left,  th* 
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Fig.  12.  Schematic  of  subsonic  curved  wall  experiment,  showing  detail  of  the 
contoured  outer  wall. 
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Fig.  13.  Profiles  of  -W/u.* 

■  .  station  1:  W,  station  3; 
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□  ,  station  9. 
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